Radiative neutralino production in low energy supersymmetric Models. II. 

The case of beam polarization 
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We study the production of the lightest neutralinos in the radiative process e + e _ — > XiX?7 m l° w 
energy supersymmetric models for the International Linear Collider energies with longitudinally 
polarized electron and positron beams. For this purpose we consider the case of nonminimal su- 
persymmetric standard model as well as the case of minimal supersymmetric standard model. At 
the first stage of a linear collider, with i/s = 500 GeV, the radiative production of the lightest 
neutralinos may be a viable channel to study supersymmetric partners of the Standard Model par- 
ticles, especially if the other supersymmetric particles are too heavy to to be pair-produced. We 
consider in detail the effect of beam polarization on the production cross section. We compare and 
contrast the dependence of the signal cross section on the parameters of the neutralino sector of the 
nonminimal and minimal supersymmetric standard model when the electron and positron beams are 
longitudinally polarized. In order to assess the feasibility of experimentally observing the radiative 
neutralino production process, we consider the background to this process coming from the Stan- 
dard Model process e + e~ — > isvj with longitudinally polarized electron and positron beams. We 
<**] , also consider the supersymmetric background to the radiative neutralino production process coming 

. from the radiative production of the scalar partners of the neutrinos (sneutrinos) e + e~ — ¥ uv*^, 

with longitudinally polarized beams. This process can be a a background to the radiative neutralino 
production when the sneutrinos decay invisibly. 

PACS numbers: 11.30.Pb, 12.60. Jv, 14.80.Ly 

O I. INTRODUCTION 

<U> '. 

Supersymmetry is at present one of the most favored ideas for physics beyond the Standard Model (SM) 0,0- 
A particularly attractive implementation of the idea of supersymmetry is the Minimal Supersymmetric Standard 
■ Model (MSSM) obtained by introducing the supersymmetric partners of the of the SM states, and introducing 
an additional Higgs doublet, with opposite hypercharge to that of the SM Higgs doublet, in order to cancel the 
gauge anomalies and generate masses for all the fermions of the Standard Model @, H| . Supersymmetry must 
\ obviously be a broken symmetry. In order for broken supersymmetry to be effective in protecting the weak 
\Q • scale against large radiative corrections, the supersymmetric partners of the SM particles should have masses 
of the order of a few hundred GeV. Their discovery is one of the main goals of present and future accelerators. 
In particular, a e + e~ linear collider with a high luminosity L = 500 fb _1 , and a center-of-mass energy of 
\fs = 500 GeV in the first stage, will be an important tool in determining the parameters of the low energy 
supersymmetric model with a high precision [5H9|. Furthermore, polarization of the electron (and positron) 
beam can enhance the capability of such a linear collider [lol ] in unravelling the structure of the underlying 
supersymmetric model. 

In the minimal supersymmetric standard model the fermionic partners of the two Higgs doublets (H\,H2) 
mix with the fermionic partners of the gauge bosons to produce four neutralino states Xi ? i — 1> 2, 3, 4, and two 
chargino states xf> 3 = 1)2. In the MSSM the lightest neutralino is favored to be the lightest supersymmetric 
particle, and assuming, i?-parity (R p ) conservation, is absolutely stable. The neutralino states of the minimal 
supersymmetric standard model with R p conservation have been studied in great detail, because the lightest 
neutralino, being the lightest supersymmetric particle (LSP), is the end product of any process involving 
supersymmetric particle in the final state. 

There are alternatives to the MSSM, an elegant one being the model with an additional chiral electroweak 
gauge singlet Higgs superfield S which couples to the two Higgs doublet superfields Hi and H 2 via a dimen- 
sionlcss trilinear term XH1H2S in the superpotcntial. This model can solve the the [i problem of the MSSM in 
a natural manner. When the scalar component of the singlet superfield S obtains a vacuum expectation value, 
a bilinear term XH1H2 < S > involving the two Higgs doublets is naturally generated. Furthermore, when 
this scalar component of the chiral singlet superfield S acquires a vacuum expectation value of the order of the 
SU(2)l x U(l)y breaking scale, it gives rise to an effective value of fi (fJ. e ff = X < S >= Xx) of the order of the 
electroweak scale. However, the inclusion of the singlet superfield leads to additional trilinear superpotential 
coupling (k/3)S 3 in the model, the so called nonminimal, or next-to-minimal [IlT4l7l ]. supersymmetric standard 
model (NMSSM). The absence of H\Hi term, and the absence of tadpole and mass couplings, S and S 2 in 
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the NMSSM is made natural by postulating a suitable discrete symmetry [1^, [l9j]. The NMSSM is attrac- 
tive on account of the simple resolution it offers to the // problem, and of the scale invariance of its classical 
action in the supersymmetric limit (20| . Since no dimensional supersymmetric parameters are present in the 
supcrpotcntial of NMSSM, it is the simplest supersymmetric extension of the Standard Model in which the 
electroweak scale originates from the supersymmetry breaking scale only. Its enlarged Higgs sector may help in 
relaxing the fine-tuning and little hierarchy problems of the MSSM, thereby opening new perspectives for the 
Higgs boson searches at high energy colliders [2l|, , and for dark matter searches [23| . In the nonminimal 
supersymmetric standard model the mixing of fermionic partners of Higgs and gauge bosons [24T-I261] produces 
five neutralino states y°, i — 1)2,3,4,5, and two chargino states xf : , j — 1,2. Furthermore, because of the 
presence of the fermionic partner of the singlet Higgs boson, the neutralino states can have an admixture of this 
SU(2)l x U(1)y singlet fcrmion, thereby affecting the phenomenology of the neutralinos in the nonminimal 
supersymmetric standard model. 

The lightest neutralino state (y? ) of MSSM or NMSSM, being typically the LSP, is stable and therefore, a 
possible dark matter candidate [27], [28| . Since the neutralinos are among the lightest particles in low energy 
supersymmetric models, they are expected to be the first states to be produced at the colliding beam experi- 
ments. At an electron-positron collider, such as the International Linear Collider (LLC), the lightest neutralino 
can be produced in pairs 

e + + e- ->x?+X?- (1-1) 

This process proceeds via Z boson and selectron exchange [29], [3(3 ■ In collider experiments the lightest neutralino 
escapes detection. In such a situation the production of the lightest neutralino pair 1|) is invisible. Therefore, 
we must look for the signature of neutralinos in the radiative process 

e~ +e + ^y? + y? + 7- (1.2) 

Despite this process being suppressed by the square of the electromagnetic coupling, it might be the first process 
where the lightest supersymmetric states could be observed at the e + e~ colliders. The signal of the radiative 
process p.2[) is a single high energy photon with the missing energy carried away by the neutralinos. The 
process (|I.2[) has been studied in detail in the minimal supersymmetric standard model [31H4l| . Some of these 
studies underline the importance of longitudinal [3{|, and even transverse beam polarizations. On the other 
hand, the signature "photon plus missing energy," that arises in the process (|I.2[) has been studied in detail by 
different LEP collaborations [42l - l46l |. Furthermore, the radiative neutrino process e + e~ — > vv^ in the SM is 
the leading process with this signature, for which the cross section depends on the number N v of light neutrino 
species [47j | . The LEP collaborations have found no deviations from the SM prediction, and, therefore, only 
bounds on the masses of supersymmetric particles have been set [42T - I441 . |46| . For a review of the experimental 
situation, see Ref. 

Most of the theoretical studies on radiative neutralino production in the literature have been carried out in 
the framework of the minimal supersymmetric standard model. This includes calculations relevant to ILC with 
a high center-of-mass energy, high luminosity and longitudinally polarized beams, as well as study of the SM 
background from the radiative neutrino production 

e+e~^^ + P + 7, (1.3) 

and the supersymmetric background from radiative sneutrino production 

e + e~~ -> v + v* + 7. (1.4) 

The discovery potential of ILC may be significantly enhanced [4!| if both beams are polarized, particularly if 
other SUSY states like heavier neutralino, chargino or even slepton pairs are too heavy to be produced at the 
first stage of the ILC at y/s = 500 GeV. 

In a previous paper (referred to as paper I) we have carried out a detailed study of the radiative process (|I.2[) in 
the nonminimal supersymmetric model and compared the predictions with those of the minimal supersymmetric 
standard model [50(. This study was carried out for unpolarizcd electron and positron beams. In this paper 
we continue this study and consider the radiative process (|I.2[) in the nonminimal supersymmetric standard 
model with polarized electron and positron beams to understand in detail if the signal can be enhanced by 
the use of polarized beams. Furthermore, the SM background photons from radiative neutrino production 
process (11.31) with beam polarizations will have to be taken into account for a proper analysis of the radiative 
neutralino production process p. 21) . Beam polarizations could enhance the signal photons for the process p.2[) 
for NMSSM and reduce those from the SM background at the same time, which could lead to the enhancement 
of the statistics. We will also consider supersymmetric background photons from radiative sneutrino production 
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process <jl.4|) with polarized beams. This is important if sneutrino production is kinematically accessible and 
if the sneutrino decay is invisible. We will compare and contrast the results obtained for NMSSM with those 
for the minimal supersymmetric standard model with polarized beams. This will include the signal for the 
radiative neutralino process, and the dependence of the cross sections on the parameters of the neutralino 
sector. This comparison will allow us to assess the feasibility of observing the radiative neutralino process for 
the most popular low energy supersymmetric models at a e + e~ collider. 

The plan of the paper is as follows. In Sec. HH we calculate the cross section for the signal process (|I.2|) 
in the nonminimal supersymmetric standard model, and compare it with the corresponding cross section in 
the minimal supersymmtric standard model, for unpolarized and polarized electron and positron beams. In 
order to calculate the cross section in the NMSSM, we fix the parameter space that we use in our calculations. 
This is done by using various theoretical and experimental constraints on the parameter space of NMSSM. In 
particular, we constrain the values of the trilinear superpotcntial parameters A and k which enter the neutralino 
mass matrix of the NMSSM. We also describe the phase space for the signal process as well as the cuts on 
outgoing photon angle and energy that we use to regularize the infrared and collinear divergences in the tree 
level cross section. We then analyze numerically the dependence of the cross section on the parameters of the 
neutralino sector, and on the selectron masses, for unpolarized and polarized beams. Here we also calculate 
the photon energy distribution for the radiative production of the second lightest neutralino in the NMSSM 
and compare it with the corresponding distribution for the lightest neutralino for unpolarized and polarized 
beams, respectively. 

In Sec. IIIII we analyze the backgrounds to the radiative neutralino process (|I.2|) with polarized beams. This 
includes the background from SM process (|I.3[) . as well as from the supersymmetric process (|I.4[) for the case 
of polarized beams. In Sec. IIVI we study in detail the beam polarization dependence of the cross section for 
the radiative neutralino production process as well as for the backgrounds processes. Here we also consider 
the statistical significance for measuring the excess of photons from radiative neutralino production over the 
backgrounds, and calculate this quantity for NMSSM, and compare it with the corresponding results in MSSM. 
We summarize our results and conclusions in Scc.|V] Our notations and results on the neutralino mass matrices 
and couplings are summarized in Appendix [A] 

II. RADIATIVE NEUTRALINO PRODUCTION 

The Feynman diagrams for the radiative process 

e-(px) + e+( P2 ) -> xUh) + xlih) + T (g), (III) 

are shown in Fig. [TJ where the symbols in the brackets denote the four momenta of the respective particles. 
In NMSSM, and in MSSM, this process proceeds at the tree level via t- and w-channcl exchange of right and 
left selectrons <?r,l, and via Z boson exchange in the s channel. In order to calculate the cross section for 
the radiative production of neutralinos we need to compute the couplings of the neutralinos to electrons, to 
the scalar partners of electrons, the selectrons, and to Z° bosons. We summarize these couplings for MSSM 
and NMSSM [5l|, [52[ in Appendix [A] As can be seen in Appendix [AJ the couplings of the lightest neutralino 
are determined by the corresponding elements of the neutralino mixing matrix (iVy or N^). For numerical 
calculation of the radiative neutralino cross section in the MSSM, we have chosen to work with the parameters in 
the Snowmass Points and Slopes (SPS la) scenario (53|. The parameters of the SPS la scenario are summarized 
in Table HI However, since in the SPS la scenario the value of the parameters /i and M 2 are fixed, we shall 
use a different set of parameters to study the dependence of the neutralino mass and the radiative neutralino 
production cross section on \i and M 2j and on the selectron masses. This set of parameters is shown in Table [TT1 
We shall call this set of parameters as the MSSM electroweak symmetry breaking scenario (EWSB) [54[. As in 
paper I, for the NMSSM we use a set of parameters that is obtained by imposing theoretical and experimental 
constraints on the parameter space of the NMSSM. The parameters that enter the neutralino mass matrix of 
the NMSSM are, apart from Mi and M 2 , tan/3, fx (= A < S >= Xx), A and k. For Mi, M 2 and M3 we use the 
values which are consistent with the usual GUT relation Mi/ai = AI 2 /a 2 = M^/a^. We note that for the MSSM 
in SPS la scenario, the value of the parameter tan/3 = 10. In order to remain close to the SPS la scenario of 
MSSM, we have chosen for our numerical calculations in NMSSM values of tan /? = 10, whereas the rest of the 
parameters are chosen in such a way that the lightest Higgs boson mass, the lightest neutralino mass and the 
lightest chargino mass satisfy the present experimental lower limits. We have also imposed on the parameter 
space of NMSSM, the theoretical constraint that there is no charge and color breaking global minimum of the 
scalar potential, and that a Landau pole does not develop below the grand unified scale (Mqut ~ 10 16 GeV). 
The consequence of imposing these constraints on the parameter space of NMSSM, and the resulting masses for 
various particles for a particular choice of input parameters is summarized in Table IIIII Since the neutralino 
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FIG. 1: Feynman diagrams for the radiative production of lightest neutralinos in the process e e 

TABLE I: Input parameters and resulting masses for various states in the MSSM SPS la scenario 
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tan P = 10 


Q = 


= 100 GeV 


m 1/2 


= 250 GeV 
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= 100 GeV 


A = 


-100 GeV 


m x o = 97 GeV 




= 180 GeV 


mg R 


= 136 GeV 




= 185 GeV 


mh = 


110 GeV 


m v o = 180 GeV 

*2 




= 379 GeV 


mg L 


= 195 GeV 


m H 


= 396 GeV 


rriA = 


395 GeV 



TABLE II: Input parameters and resulting masses of various states in MSSM EWSB scenario. 



tan ft = 10 


M = 


149 GeV 


Mi 


= 150 GeV 


M 2 


= 300 GeV 


M 3 = 


= 1050 GeV 


M A = 242 GeV 


At = 


3000 GeV 


A b -- 


= 3000 GeV 




= 1000 GeV 






m x o = 108 GeV 

m v o = -160 GeV 

x 2 


m A 


= 135 GeV 


mg R 


= 137 GeV 




= 187 GeV 


m h 


= 118 GeV 


= 328 GeV 


mg L 


= 197 GeV 


m H 


= 243 GeV 


m,A 


= 242 GeV 



mass matrix depends on the parameters A and k, it is useful to study the possible values of these parameters, 
with all other parameters fixed, which satisfy the phcnomenological and theoretical constraints discussed above. 
In Fig. [5] we show a plot of A versus k, with all other input parameters fixed as in Table IIII1 and with the 
lightest neutralino, the lightest Higgs boson, and the lightest chargino masses as in Table iLUl with a variation 
of less than 5%. Fig. [2] shows the range of A and k values that are consistent with all the constraints discussed 
above for the set of input parameters in Table IIIII We note that for the set of input values in Table IIII1 values 
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TABLE III: Input parameters and resuling masses of various states in NMSSM. 



tan P = 10 


M = 


149 GeV 


Mi 


= 150 GeV 
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300 GeV 


M 3 = 1050 GeV 


M A = 242 GeV 
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A b ~- 


= 3000 GeV 


A T = 


1000 GeV 


A = 0.54 


k = 0.45 


A A = 


880 GeV 


A K 


= 10 GeV 


rng R = 


= 137 GeV 


m SL = 197 GeV 


m x o = 94 GeV 
m v o = -160 GeV 


m x± 


= 135 GeV 


rriv e 


= 187 GeV 


m h = 


122 GeV 




m ± 


= 328 GeV 


m H2 


= 242 GeV 


m,H 3 = 


1313 GeV 






of A < 0.4, with k < 0.22, lead to an unphysical global minimum. On the other hand, values of A > 0.57, with 
k > 0.45, lead to a Landau pole below the GUT scale. Thus, the allowed values of A and k, for the given set 
of input parameters, and for the fixed masses of lightest neutralino, the lightest Higgs boson, and the lightest 
chargino, as in Table HTT1 lie in a narrow range 0.4 < A < 0.57 for 0.22 < k < 0.45. For definiteness, we have 
chosen to work with the values of A = 0.54 and K = 0.45 in this paper. These values correspond to the peak 
in the A versus K plot in Fig. [2] For the parameters of Table IIIII the composition of the lightest neutralino in 
NMSSM is given by 

N[j = (0.48, - 0.23, 0.57, - 0.55, 0.30). (II.2) 



From the composition (|II.2|) . we see that the lightest neutralino has a sizable singlet component, thereby 
changing the neutralino phenomenology in the NMSSM as compared to MSSM. For comparison, we also show 
the particle content of the lightest neutralino in MSSM 



Nij = (0.6, -0.21, 0.61, -0.47), (II.3) 



for the parameter set in Table ILT1 In Fig 
neutralino in NMSSM in the \i - M 2 p 



3(a) we have plotted the constant contour plots for the mass of lightest 
ane. We emphasize that the choice of fi and M2 values in this plot 
have been taken to be consistent with phenomenological and theoretical constraints as described above. For 
comparison, we have also plotted the corresponding contour plots for MSSM in Fig. |3(b)| with parameters as 
in Table HU 



A. Cross Section for the Signal Process 



In NMSSM, and in MSSM, the process proceeds at the tree level via t- and it-channel exchange of right 

and left selectrons Rr^l, and via Z boson exchange in the s channel. The photon is radiated off the incoming 
beams or the exchanged selectrons. The corresponding Feynman diagrams are shown in Fig.[T] The differential 
cross section for ( III.1[) can be written as 
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FIG. 3: (a) Contour plots of constant lightest neutralino mass m o in fj, -M2 plane for NMSSM; (b) for MSSM. 
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FIG. 4: (a) Photon energy distribution for the radiative neutralino production for NMSSM (red solid line), for 

MSSM EWSB (green dashed line) and MSSM SPS la (blue dashed line) at ^/s = 500 GeV with (P e - , P e+ ) =(0, 0); (b) 
with (P e -, P e+ ) =(0.8,-0.6). 



where pf and Ef denote the final three- momenta ki, k2, q and the final energies E Xl 
neutralinos and the photon, respectively. The squared matrix element \A4\ 2 in ( III.4[) can be written as [3 



E X2 , and E 1 of the 



i<3 



(II.5) 



where are squared amplitudes corresponding to the Feynman diagrams in Fig. [T] A sum over the spins of 
the outgoing neutralinos, as well as a sum over the polarizations of the outgoing photon is included in T^j. We 
have included the longitudinal beam polarizations of electrons, P e -, and positrons, P e + , with —1 < P e ± < +1, 
while calculating the cross section for the the radiative neutralino production process. The phase space in 
is described in 13311 . 



( HL4| 



Numerical Results 



We have calculated the squared amplitudes and the tree-level cross section for radiative neutralino production 
p.2j) . and the background from radiative neutrino and sneutrino production, (|I.3|) and p.4[) . with polarized 
electron and positron beams, using the program CALCHEP [54|. We note that when integrating the squared 
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FIG. 5: (a) Total energy y/s dependence of the cross sections a for radiative neutralino production e + e — > XiXil f° r 
NMSSM (red solid line) and for MSSM EWSB scenario (green dashed line)and MSSM SPS la (blue dashed line) with 
(P e -, P e +)={0, 0), (b) with (P e _, P e +) =(0.8, - 0.6). 



amplitude for the radiative neutralino production, the s — t interference terms cancel the s — u interference 
terms due to a symmetry in these channels, due to the Majorana properties of the neutralinos [39|. The tree 
level cross sections have infrared and collinear divergences, which need to be regularized [33j ]. To do this we 
define the fraction of the beam energy carried by the photon as x — E 7 / E^ cam , where y/s = 2£'boam is the 
center of mass energy, and £" 7 is the energy carried away by the photon. We then impose the following cuts 
on x, and on the scattering angle B 1 of the photon [49l ]: 



m 2 

0.02 < x < 1 



'beam 



-0.99 < cos6L < 0.99. 



(II.6) 
(II.7) 



The lower limit on x in pi. 6ft corresponds to a photon energy E~ t = 5 GeV for the center of mass energy 
y/s = 500 GeV. The upper limit of (1 — m? / E^ eam ) on x corresponds to the maximum energy that a photon 
can carry in radiative neutralino production. 

In order to implement the cuts on the photon energy in the calculation of the cross sections, we have taken 
the mass of the lightest neutralino in NMSSM to be m x a = 94 GeV for the parameter set shown in Table IIIII 
For MSSM SPSla, we take m Y o = 97 GeV and for MSSM EWSB m v o = 108 GeV. 

We note that for values of y/s = 500 GeV, and for m x o > 94 GeV, this cut reduces a substantial amount of 
the on-shcll Z boson contribution to radiative neutrino production process. 



2. Photon Energy (E 1 ) Distribution and Total Beam Energy (y/s) Dependence 

Using the procedure described above, we have calculated the energy distribution of the photons from radiative 
neutralino production in NMSSM, in MSSM SPS la, and in MSSM EWSB for both unpolarized and polarized 
electron and positron beams, respectively. These are shown in Fig. 01 where we compare the energy distribution 
of the photons in these models. In Fig. [5] we show the total beam energy y/s dependence of the cross sections for 
NMSSM, and for MSSM EWSB and MSSM SPS la, respectively. We note that the photon energy distribution 
and the total cross section for radiative neutralino production in NMSSM and in MSSM EWSB are very close 
to each other, and are smaller as compared to what one obtains in MSSM SPS la scenario. 



3. Dependence on /j, and M2 



Since the neutralino mass matrix, and hence the lightest neutralino mass, depends on [i and M2, it is 
important to study the dependence of the radiative neutralino cross section on these parameters. In the 
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FIG. 6: (a) Total cross-section a for the radiative neutralino production versus \x for NMSSM (red solid line) and for 
MSSM EWSB scenario (green dashed) at = 500 GeV with (P e _, P e+ )=(0, 0); (b) with (P e _, P e+ ) =(0.8, - 0.6). 



nonminimal supcrsymmetric standard model, fi (= A < S >= Xx) and M 2 arc independent parameters. We 
have, therefore, studied the cross section a(e + e~ — > x?Xi7) as a function of /x and M 2 independently. In Fig.[S] 
we show the fi dependence of the total cross section for the radiative production of neutralinos for NMSSM 
as well as MSSM EWSB. We recall that in the MSSM SPS la scenario these parameters are fixed. As is 
seen from Fig. \E[ the total cross section increases with /z. The plot of total cross-section versus /i in Fig. [6] is 
plotted in the range fj, € [120, 170] GeV in NMSSM and in MSSM EWSB. Note that the parameter values are 
chosen so as to avoid color and charge breaking minima, absence of Landau pole, and the phenomenological 
constraints on different particle masses. Furthermore, in Fig. [7] we show the M 2 dependence of the total cross 
section for radiative neutralino production for NMSSM and MSSM EWSB. The total cross-section decreases 
with increasing value of M 2 . The graph of total cross-section versus M 2 in Fig. [7] is plotted for the interval 
M 2 G [150,450] GeV in NMSSM and in MSSM EWSB so as to satisfy the theoretical and phenomenological 
constraints described above. From Figs. |6]and[jjwe note that the cross section is significantly enhanced when 
the electron and positron beams are polarized as compared to the case when the beams are unpolarizcd. 




160 280 400 160 280 400 

M 2 (GeV) M 2 (GeV) 



(a) (b) 

FIG. 7: (a) Total cross section a for the radiative neutralino production versus Mi for NMSSM (red solid line) and for 
MSSM EWSB scenario (green dashed) at s/s = 500 GeV with (P~ , P+) = (0, 0); (b) with (P e - , i P e +) =(0.8, - 0.6). 
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4- Dependence on selectron masses 



The cross section for radiative neutralino production er(e + e~ — > XiXil) proceeds mainly via right and left 
selectron crl exchange in the t and u-channels. In the NMSSM and MSSM EWSB, the selectron masses 
are independent parameters. In Fig. [5] and Fig. [5] we show the dependence of total cross section of radiative 
neutralino production on the left and right selectron masses. The cross section is not very sensitive to the 
selectron masses for both models. Furthermore, the total neutralino production cross section is smaller in 
NMSSM as compared to MSSM EWSB as a function of left as well as right selectron masses. Again we note 
that the cross sections, as a function of selectron masses, are larger in the case of polarized beams as compared 
to the unpolarized case. 
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FIG. 8: (a) Total cross section a for the radiative neutralino production versus mg t for NMSSM (red solid line) and 
for MSSM in the EWSB scenario (green dashed) at ^ = 500 GeV with (P e - , P e+ )=(0, 0), (b) with (P e _ , P e +) =(0.8, 
- 0.6). 
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FIG. 9: (a) Total cross section a for the radiative neutralino production versus mg R for NMSSM (red solid line) and 
for MSSM in the EWSB scenario (green dashed) at <Js = 500 GeV with (P e - , P+) = (0, 0), (b) with (P e - , P e+ ) =(0.8, 
- 0.6). 



5. Photon energy (E^) distribution for the production of the second lightest neutralino 



The cross section for the production of the lightest neutralino in NMSSM is relatively small compared with 
the corresponding cross section for the lightest neutralino in the MSSM SPSla. It may, therefore, be useful 
to consider the radiative production of the second lightest neutralino in the NMSSM. For the parameter set of 
Table UTTl the composition of the second lightest neutralino in NMSSM is given by 



NL = (0.87, 0.21, - 0.22, 0.34, - 0.19). 



(II.8) 
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125 250 125 250 



E y (GeV) E y (GeV) 



(a) (b) 



FIG. 10: (a) Photon energy distribution for the radiative production of second lightest neutralino in NMSSM (red 

solid line), and the lightest neutralino for NMSSM (green dashed line) at ,/a = 500 GeV with (P e _, P e+ ) =(0, 0); (b) 
with (P e -, P e+ ) =(0.8, - 0.6). 



Wc have calculated the photon energy distribution for the radiative production of the second lightest neutralino 
in NMSSM for the set of parameters shown in Table Hill This is shown in Fig. [TO] For comparison we have also 
shown the photon energy distribution for the radiative production of the lightest neutralino in NMSSM. We 
see that the cross section for the production of the second lightest neutralino is much smaller than the cross 
section for the lightest neutralino at photon energy greater than 140 GeV. However, at lower photon energies 
the photon energy distribution for the second lightest neutralino is significantly larger, both for unpolarized as 
well as polarized beams. 
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FIG. 11: Feynman diagrams contributing to the radiative neutrino process e + e~ — > v&j. 
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125 250 125 250 



E y (GeV) E y (GeV) 

(a) (b) 

FIG. 12: (a) The photon energy distribution for the radiative neutrino process e + e~ -> vi>j at yfs — 500 GeV 
with (P e -, P+)=(0, 0), (b) with (P e _, P e+ ) =(0.8, - 0.6). 



III. BACKGROUND PROCESSES 
A. The Neutrino Background 

The major background to the radiative neutralino production ( comes from the SM radiative neutrino 
production process [i(J 53, l56l - [58j 

e+ + e~ -> v t + D e +7 , i = e,n,T. (HIT) 

In this process ^ e are produced via ^-channel W boson exchange, and v e ,^,r via s-channel Z boson exchange. 
We Feynman diagrams contributing to the process are shown in Fig. [TTJ 

The background photon energy distribution j^- and y/s dependence of the cross section a for radiative 

neutrino production e + e~ — > with polarized electron and positron beams is the same for both NMSSM 
and MSSM. As shown in Fig.[T5]the photon energy distribution from the radiative neutrino production peaks at 
E 1 = (s — m z ) I '(2y/s) ss 242 GeV because of the radiative Z production( v / s > mz). This photon background 
from radiative neutrino production can be reduced by imposing an upper cut on the photon energy x max = 
E™ ayL /E bcam = 1 - m 2 x0 /El cam GeV in NMSSM, see Eq. (ITLoTl . which is the kinematics] endpoint £™ ax 215 
GeV of the energy distribution of the photon from radiative neutralino production 

m 2 x o =\(a- 2^~sE^) . (III.2) 

In order to achieve this, one would have to separate the signal and background processes. This would be possible 
if the neutralino is heavy enough, such that the endpoint is removed from the Z° peak of the background 
distribution. See also [49||. In Fig. [T3]we show the -^/s dependence of the total radiative neutrino cross section 
with unpolarized and polarized electron and positron beams without imposing upper cut on the photon energy. 
In Section HVl the upper cut on the photon energy i?™ ax = 214.7 GeV is used for the calculation of cross section 
for the radiative neutrino production. 



B. The Supersymmetric Background 

Apart from the SM background coming from (|III.1[) . the radiative neutralino production (|II.1[) has a back- 
ground coming from the supersymmetric sneutrino production process [10, H^] 

e+ + e" -> v t + &t +7, l = e,n,r. (III.3) 

The lowest order Feynman diagrams contributing to the process (|III.3[> are shown in Fig. [TU This process 
receives i-channel contributions via virtual charginos for z/ e i/*-production, as well as s-channel contributions 
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FIG. 13: (a) Total energy yfs dependence of the cross sections a for radiative neutrino cross section a(e + e — > wy) 
with (P e _, P e+ )=(0, 0), (b) with (P e _, P e+ ) =(0.8, - 0.6). 

from Z boson exchange for v e ^^ r v* ^ r -production. In Fig. 1151 we show the photon energy distribution -j^- for 

radiative sneutrino production e + e~ — > vv*~{ at ^fs = 500 GeV with unpolarized and longitudinally polarized 
electron and positron beams. The corresponding total cross section for the radiative sneutrino production is 
shown in Fig. 1161 
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FIG. 14: Feynman diagrams contributing to the radiative sneutrino production process e + e 



vv 7. 



Radiative sneutrino production ( IIII.3|) can be a major supersymmetric background to neutralino production 
( III.1[) if sneutrinos decay mainly invisibly, e.g. via v — > xfv. This leads to so called "virtual LSP" scenario (ioj . 



However, if kinematically allowed, other visible decay channels like v — > reduce the background rate 

from radiative sneutrino production. For example in the SPS la scenario [HD, H3| of the MSSM we have 
BR(£> e -> x\u e ) = 85%. 

Furthermore, neutralino production e + e~ — > X1X2 followed by subsequent radiative neutralino decay [6l[ 
X2 ~> Xil is also a potential background. However, significant branching ratios BR(x§ —> Xil) > 10% are 
obtained only for small values of tan/3 < 5 and/or Mi ~ M2 [HI, [62|, |H| . Thus, we neglect this background, 
detailed discussions of which can be found in Refs. [62| - [64{ . 
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FIG. 15: (a)The photon energy distribution -4§- for the radiative sneutrino production e + e 
GeV with (P e _, P e +)=(0, 0), (b) with (P e _, P e+ ) =(0.8, - 0.6). 
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FIG. 16: (a)Total energy y/s dependence of the radiative sneutrino production cross section a(e + e — > £z/*7)with 
(P e _, P e+ )=(0, 0), (b) with (P e _, P e+ ) =(0.8, - 0.6). 



IV. THE EFFECT OF LONGITUDINAL BEAM POLARISATION 



In this Section we study in detail the effect of longitudinal beam polarization on the cross section for 
radiative neutralino production in NMSSM, and MSSM, in electron-positron collisions. We shall also study the 
beam polarization dependence of the background processes. The cross section with polarized e^- beams, with 
polarization P e ± (\P e ± | < 1), can be written as 



a(P e -,P e+ ) = 2J2 



1 + P P -X\ /l-P„+\ 



(IV.l) 



where a\ is the helicity dependent cross section, and where A is the helicity (A = ±1) of the initial particles (e^). 
The signal radiative neutralino production cross section can be enhanced, and the background can be reduced 
by using a positively-polarized e~ beam (P e - > 0) and a negatively-polarized e + beam (P e + < 0). 

In the NMSSM, and MSSM, the radiative neutralino production process proceeds proceeds mainly via the 
exchange of right selectrons en- This is because, for the parameter choices that we use, the neutralino has a 
significant bino component in these models, and the coupling to the right selectron is significantly larger as 
compared to the left selectrons e^. This leads to the contribution from the right selectron exchange to the 
cross section which is an order of magnitude larger as compared to the left selectron exchange. Furthermore, 
compared to the right selectron exchange, the the contribution from left selectron exchange is suppressed due 
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TABLE IV: Cross sections er, significance S, and signal to background ratio r for different beam polarizations (P e - \P e + ) 
for MSSM SPS la at -s/s = 350 GeV, with tan/3 = 10 and £ = 500 fb -1 . 





(0]0) 


(0.8|0) (0.8|-0.3) (0.8|-0.6) (0.9|0) (0. 


9| -0.3) (0.9| -0.6) 




20 fb 


35 fb 


46 fb 


56.4 fb 


37 fb 


18.4 fb 


60 fb 


cr(e + e~ — > vv^j) 


2991 fb 


793 fb 


685 fb 


579 fb 


518 fb 


501 fb 


484 fb 


S 


8 


27 


38 


50 


35 


46 


57.5 


r 


0.7% 


4.4% 


6.7% 


9.7% 


7.1% 


9.6% 


12.4% 



TABLE V: Cross sections a, significance S, and signal to background ratio r for different beam polarizations (P e - \P e +) 
for MSSM SPS la at y/s = 500 GeV, with tan /3 = 10 and C = 500 fb" 1 . 



(P e -\Pe+) 


(0]0) 


(0.8|0) (0 


8| -0.3) (0 


8| — 0.6) (0.9|0) (0.91-0.3) (0.9] - 0.6) 


a( e + e- -> XiXh) 


23 fb 


40 fb 


52 fb 


64 fb 


42 fb 


55 fb 


67 fb 


a(e + e~ —> vv^f) 


2942 fb 


597 fb 


423 fb 


250 fb 


304 fb 


218 fb 


133 fb 


S 


9.4 


35.4 


53.3 


80.8 


50.5 


74.4 


106 


r 


0.8% 


6.7% 


12.3% 


25.6% 


13.8% 


25.2% 


50.4% 



TABLE VI: Cross sections a, significance S, and signal to background ratio r for different beam polarizations [P e - \P e + ) 
for MSSM SPS la at y/s = 650 GeV, with tan/3 = 10 and C = 500 fb" 1 . 



(Pe-\P e +) 


(010) 


(0.8|0) (0.8|-0.3) (0.8]-0.6) (0.9|0) (0.9) - 0.3) (0.9| - 0.6) 


o{e + e- -> XiXi7) 


18 fb 


32 fb 


42 fb 


52 fb 


34 fb 


44 fb 


55 fb 


a(e + e~ — > vv"y) 


2869 fb 


578 fb 


408 fb 


237 fb 


292 fb 


207 fb 


123 fb 


S 


7.5 


29 


44 


68 


42 


62 


92 


r 


0.63% 


5.5% 


10.3% 


22% 


11.6% 


21.3% 


44.7% 



to the fact that m,g R < m^ L . We also note that the Z boson exchange contribution to the neutralino production 
process is negligible in these models. The SM background radiative neutrino process proceeds mainly via the 
exchange of W bosons. This means that positive electron beam polarization P e - and negative positron beam 
polarization P e + will enhance the signal cross section, and at the same time reduce the background [3!|. We 
note that \P e - | > 0.8 and \P e +\ > 0.5 are designed at the International Linear Collider (6f|. 

A quantitative measure of the excess of photons from the radiative neutralino production over the SM 
background photons is the theoretical significance defined as 

S = - Ns = - as VC, (IV.2) 

y/(N 8 + N B ) V^S + CTB 

where Ns and Nb define the number of signal and background events, respectively, as and ae are the respective 
cross sections, and C is the integrated luminosity. If the theoretical significance has a value S = 1 for a signal, 
then that signal can be measured at a 68% confidence level. In addition, we must also consider the signal to 
background ratio 

r = — . (IV.3) 
aB 

For a signal to be detectable at the International Linear Collider we must have S > 1 and r > 1% . These 
estimates are expected to be rough estimates which will enable us to decide whether an excess of signal photons 
can be measured over the background photons. A detailed Monte Carlo analysis is beyond the scope of the 
present paper. 

We have studied the radiative production of lightest neutralino for three different electron-positron center 
of mass energies, namely y/s = 350 GeV, 500 GeV, and 650 GeV, respectively with longitudinally polarized 
beams for the three models— MSSM SPSla, MSSM EWSB and NMSSM using CALCHEP. For MSSM SPS 
la, we have calculated the beam polarization dependence of the signal cross section a(e + e~ —> XiXil) an d 
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the background cross section er(e + e~ — > vP~f), the significance S and signal to background ratio r at y/s = 350 
GeV, 500 GeV and 650 GeV for the input parameters as in Table HI We report the values of cross-sections 
er(e + e - — > XiXil) an d cr(e+e~ — > vv'f), S and r for a specific set of beam of polarizations (P e - |P e +) = (0|0), 
(0.8|0), (0.8|-0.3), (0.8|-0.6), (0.9|0), (0.9| — 0.3), and (0.9| —0.6) in Tables IT71 El and EH respectively. We 
observe that there is a large enhancement in the value of S = 106 and r = 50.4% for (P e - \P e +) = (0.9 1 — 0.6) 
compared to all other chosen set of polarization values for y/s = 500 GeV. For a previous study of the beam 
polarization effects for the radiative neutralino production in MSSM, see Ref. [6a |. 

Similarly, in the Tables IVTI1 IVIII1 and IIX1 we discuss the values of cross-sections a(e + e~ —> XiXil) an d 
<r(e + e~ — > uP^f), S and r for a specific set of beam polarization values (P e - \P e +) = (0|0), (0.8|0), (0.8 1 — 0.3), 
(0.8| - 0.6), (0.9|0), (0.9| - 0.3), and (0.9| - 0.6) for the case of MSSM EWSB. There is again a significant 
increase for the value of S = 12.5 and r = 5% for (P e - \P E +) = (0.9| —0.6) compared to other chosen polarization 
values for y/s = 500 GeV. However, the enhancement is relatively samll as compared to the MSSM SPS la 
scenario. 

Finally, in Tables [Xj IXI1 IXII1 we show the results for the case of NMSSM. We get the similar pattern but 
the value of S — 7 and r = 2.8% for (P e - \P e +) = (0.9 1 — 0.6) is considerably smaller as compared to that of 
S = 106 and r = 50.4% in MSSM SPS la at y/s = 500 GeV. Thus, if the radiative production of neutralinos 
is observed at a linear collider with polarized beams, then on the basis of the observed event rate it may be 
possible to distinguish between MSSM and NMSSM as the underlying low energy supersymmetric model. 

V. SUMMARY AND CONCLUSIONS 

The nonminimal supersymmetric standard model (NMSSM) solves /i problem of MSSM in an elegant manner 
and is, thus, an attractive alternative to the MSSM. We have carried out a detailed study of the radiative 
neutralino production process e + e~ — > XiXil m NMSSM for ILC energies, and compared the results with the 
corresponding results in the MSSM for both unpolarized and polarized and e + beams. This process has 
a signature of a high energy photon and missing energy. We have obtained a typical set of parameter values 
for the NMSSM by imposing theoretical and experimental constraints on the parameter space of NMSSM. For 
the set of parameter values that we obtain in this manner, the lightest neutralino in NMSSM has a significant 
admixture of the fermionic component of the singlet chiral supcrficld S. Using this parameter set, we have 
studied in detail the radiative neutralino production cross section in NMSSM for the ILC energies for both 
unpolarized and polarized e~ and e + beams. For comparison with MSSM, we have used the MSSM SPS la 
and MSSM EWSB models. The background to this process comes from the SM process e + e~ — > vvj, as well 
as the background from the supersymmetric process e + e~ — > vv*^. All these processes have a signature of 
a highly energetic photon with missing energy. The photon energy distribution da/dE^, and the total cross 
section as a function of the total energy have been calculated for the NMSSM and for MSSM SPS la scenario 
at y/s = 500 GeV using CALCHEP package. Because of the admixture of a singlet in the lightest neutralino, 
the cross section as a function of energy for the radiative neutralino production is lower in NMSSM than in 
MSSM. We have also studied the dependence of the cross section for radiative neutralino production on the 
SU(2)l gaugino mass parameter M 2 and the Higgs(ino) mass parameter /x, as well as its dependence on the 
selectron (e^e^) masses in NMSSM, and compared it with the corresponding results in MSSM. In order to 
quantify whether an excess of signal photons, iVg, can be measured over the background photons, Nb, from 
radiative neutrino production, we have analyzed the theoretical statistical significance S = Ng / y/Ng + Nb ■ At 
the ILC, electron and positron beam polarizations can be used to significantly enhance the signal and suppress 
the background simultaneously. We have shown that the significance can then be increased almost by an order 
of magnitude, e.g. with (P e -,P e +) = (0.8,-0.6) compared to (P e -,P e +) = (0,0). A possible feed-back of 
ILC500 results could motivate the immediate use of the low-luminosity option of the ILC at yfs = 650 GeV in 



TABLE VII: Cross sections a, significance 5", and signal to background ratio r for different beam polarizations (P e - ]P C + ) 
for MSSM EWSB at^s^ 350 GeV, with tan/? = 10 and £ = 500 fb" 1 . 



(P e -lp e +) 


(0|0) 


(0.8|0) (0.8|-0.3) (0.8|-0.6) (0.9|0) (0.9| - 0.3) (0.9| - 0.6) 


°{e + e~ XiXi7) 


1.7 fb 


2.9 fb 


3.8 fb 


4.7 fb 


3.1 fb 


4 fb 


4.9 fb 


a(e + e~ — > vv^j) 


2991 fb 


793 fb 


685 fb 


579 fb 


518 fb 


501 fb 


484 fb 


S 


0.7 


2.3 


3.2 


4.3 


3 


4 


5 


r 


0.06% 


0.4% 


0.56% 


0.81% 


0.6% 


0.8% 


1.02% 
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TABLE VIII: Cross sections a, significance S, and signal to background ratio r for different beam polarizations 
(Pe- \P e +) for MSSM EWSB at Vs = 500 GeV, with tan p = 10 and C = 500 fb _1 . 



(P e -\p e+ ) 


(0]0) 


(0.8|0) (0.8|-0.3) (0.8) -0.6) (0.9|0) (0 


9| -0.3) (0.9| -0.6) 


o{e + e- -> XiXh) 


2.2 fb 


3.9 fb 


5.1 fb 


6.3 fb 


4.1 fb 


5.4 fb 


6.6 fb 


cr(e + e~ — > vv^j) 


2942 fb 


597 fb 


423 fb 


250 fb 


804 fb 


218 fb 


133 fb 


S 


0.9 


3.6 


5.5 


8.8 


5.2 


8 


12.5 


r 


0.075% 


0.65% 


1.2% 


2.5% 


1.3% 


2.5% 


5% 



TABLE IX: Cross sections a, significance 5*, and signal to background ratio r for different beam polarizations (P e - \P e + ) 
for MSSM EWSB at s/s = 650 GeV, with tan /3 = 10 and £ = 500 fb -1 . 



(Pe-\P e +) 


(0]0) 


(0.8|0) (0.8|-0.3) (0.8]-0.6) (0.9|0) (0.91-0.3) (0.9| - 0.6) 


°{e + e- -> XiXi7) 


1.9 fb 


3.4 fb 


4.4 fb 


5.4 fb 


3.5 fb 


4.6 fb 


5.7 fb 


a(e + e~ — > ui"y) 


2869 fb 


578 fb 


408 fb 


237 fb 


292 fb 


207 fb 


123 fb 


S 


0.8 


3 


4.8 


7.7 


4.5 


7 


11 


r 


0.06% 


0.6% 


1% 


2.3% 


1.2% 


2.2% 


4.6% 



TABLE X: Cross sections a, significance S, and signal to background ratio r for different beam polarizations (P e - \P e +) 
for NMSSM at s/s = 350 GeV, with tan/3 = 10 and C = 500 fb" 1 . 



(Pe-\P e+ ) 


(0]0) 


(0.8|0) (0.8|-0.3) (0.8|-0.6) (0.9|0) (0.91-0.3) (0.9| - 0.6) 




1.1 fb 


2 fb 


2.6 fb 


3.2 fb 


2.1 fb 


2.7 fb 


3.3 fb 


cr(e + e~ — > vv^j) 


2991 fb 


793 fb 


685 fb 


579 fb 


518 fb 


501 fb 


484 fb 


S 


0.45 


1.6 


2.2 


3 


2 


2.7 


3.3 


r 


0.04% 


0.25% 


0.4% 


0.55% 


0.4% 


0.54% 


0.68% 



TABLE XI: Cross sections a, significance S, and signal to background ratio r for different beam polarizations (P e - \P e +) 
for NMSSM at s/s = 500 GeV, with tan/3 = 10 and C = 500 fb" 1 . 



{P e -\p e +) 


(010) 


(0.8|0) (0.8|-0.3) (0.8)-0.6) (0.9|0) (0.91-0.3) (0.9| - 0.6) 


a(e + e~ -> Xixh) 


1.2 fb 


2.2 fb 


2.8 fb 


3.5 fb 


2.3 fb 


3 fb 


3.7 fb 


a(e + e~ — > vv^j) 


2942 fb 


597 fb 


423 fb 


250 fb 


804 fb 


218 fb 


133 fb 


S 


0.5 


2 


3 


5 


3 


4.5 


7 


r 


0.04% 


0.4% 


0.66% 


1.4% 


0.76% 


1.4% 


2.8% 



TABLE XII: Cross sections a, significance 5", and signal to background ratio r for different beam polarizations (P e - \P e + ) 
for NMSSM at yfs = 650 GeV, with tan P = 10 and C = 500 fb" 1 . 



(P e -\p e+ ) 


(0)0) 


(0.8|0) (0.8|-0.3) (0.81-0.6) (0.9|0) (0.9) - 0.3) (0.9| - 0.6) 


o{e + e~ XiXi7) 


1 fb 


1.8 fb 


2.3 fb 


2.8 fb 


1.9 fb 


2.4 fb 


3 fb 


cr(e + e" — > vv~j) 


2869 fb 


578 fb 


408 fb 


237 fb 


292 fb 


207 fb 


123 fb 


S 


0.42 


1.7 


2.5 


4 


2.5 


3.7 


6 


r 


0.03% 


0.3% 


0.56% 


1.2% 


0.65% 


1.16% 


2.44% 



order to resolve model ambiguities between NMSSM and MSSM even at an early stage of the experiment and 
outline future search strategies at the upgraded ILC at 1 TeV. In our scenarios, the signal cross section for 
(P e - \P e + ) = (0.8 1 — 0.6) is larger than 3.5 fb, the significance S > 5, and the signal to background ratio is about 
r > 1%. The background cross section can be reduced to 250 fb. Increasing the positron beam polarization 
to P e + = —0.6, both the signal cross section and the significance increase significantly. Thus the electron and 
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positron beam polarization at the ILC are essential tools to observe radiative neutralino production. 
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Appendix A: Superpotential, neutralino mass matrix and couplings 



For completeness we summarize here the couplings of the neutralinos to electrons and the scalar partners 
of electrons, the selectrons, in MSSM and in NMSSM. These couplings can be obtained from the neutralino 
mixing matrix. To obtain the neutralino mixing matrix for the MSSM, we recall that the neutralino mass 
matrix obtains contributions from part of the MSSM superpotential 

Wmssm = fiHiHa, (A.l) 

where Hi and H 2 are the two Higgs doublet chiral supcrficlds, and fi is the supersymmetric Higgs(ino) parame- 
ter. In addition to the contribution from the superpotential, the neutralino mass matrix receives contributions 
from the interactions between gauge and matter multiplets, as well as contributions from the soft supersym- 
metry breaking masses for the gauginos. Including all these contributions, the neutralino mass matrix, in the 
bino, wino, higgsino basis (— iX', — iA 3 , ipjj , ipg ) can be written as (HlL l52| 

/ Mi — mz sin 9 W cos (3 mz sin 9 W sin /3\ 

M 2 mz cos 6 W cos /3 — mz cos 6 W sin /? I . 

MSSM ~ \-m z sin 6 W cos f3 m z cos 9 W cos p -p I ' 1 ' 

\ mzshi8 w su\(3 — mz cos 9 W sin /3 —fi / 

where Mi and M 2 are the U(1)y and the SU(2)l soft supcrsymmetry breaking gaugino mass parameters, 
respectively, and tan j3 = v 2 jv\ is the ratio of the vacuum expectation values of the neutral components of 
the two Higgs doublet fields Hi and H 2 , respectively. Furthermore, mz is the Z boson mass, and 9 W is the 
weak mixing angle. We shall assume that all the parameters in the matrix M are real, in which case M is 
a real symmetric matrix and can be diagonalised by an orthogonal matrix. If one of the eigenvalues of M is 
negative, one can diagonalize this matrix using a unitary matrix N, the neutralino mixing matrix, to get a 
positive diagonal matrix [5l| : 

^Mmssm^" 1 = diag(m x o, to % o, m x o, m x o) . (A.3) 

where m x q (i = 1, 2, 3, 4) are neutralino masses arranged in order of increasing value. 
For the NMSSM, the relevant part of the superpotential is 

V^nmssm = XSH1H2 - ~S 3 , (A.4) 

where S is the Higgs singlet chiral supcrfield. In the basis (— iX', — iA 3 , ipjj , f/'lf ,ips), the neutralino mass 
matrix for the NMSSM can then be written as [H [H 



NMSSM 



/ 


Mi 





—mz sin 9 W 


cos j3 


mz sin 9 W 


sin/3 


° \ 







M 2 


m z cos 9 W 


cos/3 


—mz cos 9 W 


sin j3 







-mz sin 9 W cos (3 


mz cos 9 W cos j3 







—Xx 




-Xv 2 




mz sin W sin j5 


—mz cos 9 W sin/3 


—Xx 









—Xvi 


V 








-Xv 2 




—Xvi 




2kx / 



(A.5) 



where < S >= x is the vacuum expectation value of the singlet Higgs field. As in the case of MSSM, 
we can use a unitary matrix TV' to get a positive semidefinite diagonal matrix with the neutralino masses 
m x o (i = 1,2,3,4,5) |2i|2l: 

N'* M^mssmN'- 1 = diag(m x o, m x o, m x o, m x o m x o) . (A.6) 
The Lagrangian for the interaction of neutralinos, electrons, selectrons and Z bosons for MSSM is given by [5l| 

C = (-^- N n)fe p LXi^ + r- e - (Ni 2 + t & u9 w Nn)f e P R xle L 
cos 9 W V2 sin 9 W 



4 sin 6,„ cos 9„ 



(\Ni 3 \ 2 \Ni 4 \ 2 ) Z^l'xi 



+eZje'f [ . 1 (\ - sin 2 9 W ) P L - tan 9 W P R ] f e + h.c, (A.7) 

sin uj cos u yj v £ j 

with the electron, selectron, neutralino and Z boson fields denoted by f e , cl,r, Xii an£ i respectively, and 
Pr,l — \ (l ± 7 5 ) ■ The corresponding interaction Lagrangian for NMSSM is obtained from ( IA.7|) by replacing 
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TABLE XIII: Vertices corresponding to various terms in the interaction Lagrangian ( lA.Tf) for MSSM. In addition we 
have also shown the vertices for selectron-photon and electron-photon interactions. The vertices for the NMSSM are 
obtained by replacing Nij with iVy. 



Vertex 



Vertex Factor 



right selectron - electron - neutralino 
left selectron - electron - neutralino 
neutralino - Z° - neutralino 
electron - Z - electron 
selectron - photon - selectron 
electron - photon - electron 



•^2 iy* p 



J* (N12 + tan e w Nn)P R 

V 2 sin G w v ' 

(|AT 13 | 2 - |A 14 | 2 ) 7 M 7 5 



4 sin 6 W cos ( 



1 - sin 2 9 W ) P L - tan W P R ] 



_ sin m cos Q -in V 2 



ie(pi +P2Y 



with iVy. The different vertices following from ( IA.7[) arc shown in Table IXTlTI The couplings of the lightest 
neutralino to electrons, selectrons and Z boson are determined by the corresponding elements of the neutralino 
mixing matrix (Nij or N^). 
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